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Extent of force indeterminacy in packings of frictional rigid disks 
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Static packings of frictional rigid particles are investigated by means of discrete element simu- 
lations. We explore the ensemble of allowed force realizations in the space of contact forces for a 
given packing structure. We estimate the extent of force indeterminacy with different methods. The 
indeterminacy exhibits a nonmonotonic dependence on the interparticle friction coefficient. We ver- 
ify directly that larger force-indeterminacy is accompanied by a more robust behavior against local 
perturbations. We also investigate the local indeterminacy of individual contact forces. The prob- 
ability distribution of local indeterminacy changes its shape depending on friction. We find that 
local indeterminacy tends to be larger on force chains for intermediate friction. This correlation 
disappears in the large friction limit. 
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In packings of relatively rigid particles, elastic defor- 
mations of the grains are typically several orders of mag- 
nitude smaller than the grain size. Since this separation 
of length scales occurs it is a natural idea to investigate 
the limit case of infinite stiffness of the grains. 

It is known that jammed packings of perfectly rigid 
particles with finite friction coefficient are ^^hyperstatic^^ 
[l|, |2| • The number of equations of mechanical balance 
is smaller than the number of unknowns (components of 
the interparticle forces). This makes the problem unde- 
termined in the sense that there are many solutions that 
satisfy the equilibrium equations. Even taking constraint 
conditions, like Coulomb's limit of friction and unilater- 
ality of the contacts, into account does not help to elim- 
inate the indeterminacy of the contact forces; Thus for a 
given packing geometry the solutions define an ensemble 
of admissible force networks 5 0, |j| . iS is a convex set 
Q in the force space J-, spanned by the components of 
contact forces, and its boundaries are delimited by con- 
straint conditions. 

The ensemble has received considerable attention since 
many macroscopic properties of granular packings can be 
derived from ensemble averaging over all allowed force 
states supposing a uniform measure on S [y, 0, [a, 11, [lOj 
lll| . Furthermore, with this technique one can disentan- 
gle the effect of forces and texture of the packing. Math- 
ematically, the problem of finding the solutions of a set 
of undetermined equations and constraints is of rather 
broad interest, e.g. in metabolic networks [12.ll3j. 

The extent of force indeterminacy in 2D random pack- 
ings of perfectly rigid disks was investigated theoretically 
and numerically in [14| . The indeterminacy of each com- 
ponent of the contact forces was obtained, suggesting 
that highly undetermined contacts are located on main 
force chains. Force indeterminacy in such packings was 
also measured in [5| where it turned out that the inde- 
terminacy depends nonmonotonically on the interparti- 
cle friction coefficient due to the competition between two 



coexisting effects, the opening of the Coulomb cone angle 
and the lowering of connectivity. In Ref. [15| similar non- 
monotonic friction dependence is obtained for mechanical 
response of the granular packings to local perturbations. 
In this paper we examine whether the nonmonotonic 
friction dependence of force indeterminacy remains valid 
also when other methods are used to quantify the "size" 
of the solution set S. We measure numerically the ex- 
tent of force-indeterminacy and the mechanical response 
to local perturbations in the same packings and examine 
the relation between them. The local force indetermi- 
nacy is also studied in this work. First, we investigate 
its probability distribution, then, we compare its spatial 
pattern to that of the force chains in the packing. 

Sampling Procedure - The systems we investigate are 
2D random packings of 400 perfectly rigid disks. Peri- 
odic boundary conditions are applied in both directions, 
disk radii are uniformly distributed between 0.5 and 1, 
gravity is set to zero and the unit of the length is set to 
the maximum grain radius. Our numerical simulations 
consist of two steps which are performed with the help 
of contact dynamics algorithm [1^, [13, [lH . First we con- 
struct static configurations of particles. The initial dilute 
systems are compressed by imposing a homogeneous con- 
fining pressure Pq to get the final static packings. The 
full description of our method of constructing the homo- 
geneous packings can be found in [19.] . 

Then we explore the force ensemble: we collect force 
networks that provide static solutions for the given con- 
tact geometry and boundary conditions. We use a ran- 
dom walk method in the force space Q, [2^ starting with 
the original force network. We perturb the original force 
state and jump to a new force state in the force space J-'. 
The technique is to add random values that are chosen 
uniformly from the interval [— (i^„), (i^„)] to all compo- 
nents of the contact forces. (Fn) is the mean normal 
force calculated over the current values of contact forces. 
The perturbed force network is given as the input for the 




FIG. 1: Force indeterminacy ri, quantified with different 
metiiods ttiat are expiained in tlie text, in terms of ttie fric- 
tion coefficient fi. The inset dispiays the same plot in iog-iog 
scaie. 



Gauss-Seidel-like iterative solver of the contact dynamics 
method which lets the forces relax into a new consistent 
state. The jump is accepted if the new state is an equi- 
librium state, otherwise it is rejected. The perturbation 
and relaxation are repeated many times, always starting 
from the last equilibrium force network. In this way we 
collect 1000 admissible force networks for a given static 
packing. In order to study systematically the influence 
of the interparticle friction coefficient on the extent of 
the force indeterminacy, the constructing and sampling 
procedures are repeated for various values of the friction 
coefficient. 

Numerical Results - Next we quantify the extent of 
the force indeterminacy rj for a given packing geometry 
based on the sampled force networks. We compare here 
three different methods. Let us denote the center of the 
samples in the force space !F by {Gc} which is a force 
network with contact force vectors Gc given by 



where the average (■ ■ 



c=l,...,A^c 
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states i^ taken over all realizations 
of the force states and Nc is the number of contacts. 
One possibility to quantify the force indeterminacy is to 
measure the force fluctuations SFc around the mean force 
vector Gc at each contact c [5] : 



SFc^{{Fc-Gc 



(2) 



The force indeterminacy 771 of the whole packing is given 
by the relative fluctuation: 
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FIG. 2: The average critical force {F„it) scaled by the average 
normal contact force (Fcont) (left) and the decay exponent a 
(right) as functions of the force indeterminacy 772 for packings 
that are constructed with different friction coefficients. 



The extent of the indeterminacy could be also esti- 
mated by the Euclidean distance between randomly cho- 
sen pairs of force states p, ll3| . The probability distribu- 
tion of the distances becomes sharply peaked if iS is a high 
dimensional object. The global indeterminacy according 
to this method is defined via 



'{i{Fc}-{Fc}r), 

{GcV 



(4) 



{Fc} and {Fc}' are two different force states and (• • •)pairs 
means the average over all pairs of force states. The 
square of a force state {Fc}^ is given by 2. Fc- 

, , c 

As an alternative method [ij] , the extremal points of S 
along each axis of the force space J^ provide the following 
measure of the indeterminacy: 



(F" 
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where (• • ■)(.o,jt denotes the average over all contacts. 



Here, F"^^'^ and i<"™™ are the maximum and minimum 
values of a contact force component (either normal or 
tangential). The average (•••)^qj^ is taken over all 2 xA^c 
components of contact forces. We note that the first two 
methods [Eqs. ([3]) and ^] depend on the probability 
measure that is realized by the sampling. This is not 
the case for 773 which has a pure geometrical definition 
(provided the sampling explores the solution set). The 
question whether the sampling is uniform or not has no 
effect on the value of 773 . 

In Figure [1] we compare the values of 77 obtained by the 
three methods which, up to a constant factor, provide ba- 
sically the same behavior in the whole range of friction 
(jii K, rj2 ~ 0.15773). The nonmonotonic friction depen- 
dence, reported in [5|, is reproduced here independently 
of the quantifying method. 

Next we investigate the effect of rj on the mechanical 
response of granular packings. In Refs. [lalSlI local per- 
turbations were used to break the equilibrium structure 
of the homogeneous packings and induce motion of the 
grains. It turned out that the displacements of the par- 
ticles due to local perturbations decay as a power law of 
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FIG. 3: The possible values of the normal contact forces F„ 
(denoted by dots and mostly merged to intervals), scaled by 
the normal component of the original contact force F°^^^, are 
shown for each contact of the packing. Each plot corresponds 
to the packing with the given friction coefficient. 
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FIG. 5: The position of the contacts with original normal 
force Fn larger than 2{Fn) (0) and the position of the contacts 
with force indeterminacy A_F„ larger than {AF„) (•) in the 
packing with /^ = 0.1 (a) and /x = 10 (b). 



the distance from the perturbation point. The numerical 
experiment was repeated for several packings constructed 
with different fj, revealing that the decay exponent a and 
the critical force .Fcrit, i-e. the force needed to break the 
mechanical equilibrium, exhibit a nonmonotonic depen- 
dence on the friction with extrema at /x = 0.1 similarly 
to the behavior of rj. This similarity suggest the picture 
that a packing with larger force-indeterminacy becomes 
more stable against perturbations. Here, we test directly 
whether such a relation exists: Together with the force 
indeterminacy we determine also the response quantities 
-Fcrit and a for the same packing configurations. Since 
the different methods we used to quantify 77 are basically 
equivalent, we plot the response quantities in terms of 77, 
in Fig. [21 The same series of packings are plotted here 




FIG. 4: The probability distribution of the interval of pos- 
sible normal contact forces for several friction coefficients /i. 
The insets display semilogarithmic probability distributions 
for three different p. 



as in Fig. [TJ The lines are connecting the data points 
in the order of increasing friction. Both i^crit and a are 
strongly related to the extent of force indeterminacy, al- 
though they are not a unique function of rj. Still, very dif- 
ferent packings (with different density, connectivity and 
frictional properties) exhibit similar response properties 
if their 77 values are close to each other. 

Next we study the local indeterminacy at the level of 
individual contacts. In Fig. [3] we show the values of nor- 
mal contact forces Fn for every contacts obtained by the 
sampled realizations of force states. At each contact the 
possible values of Fn form an interval due to the convex- 
ity of the solution set S. The length of the interval AF„ 
can be estimated with help of the extrema of Fn that 
were provided by the sampling procedure. The values 
of AFn are very small in the nearly frictionless packing 
(/i = 10"®). The intervals become wider with increasing 
the friction, but only up to /x = 0.1. Beyond this point 
the local indeterminacy starts decreasing. 

Fig. [5] shows the probability distribution of AF„ for 
different friction coefficients. We find that P(AF„) is 
a monotonically decreasing function for small and large 
friction limits, but becomes broader and displays a peak 
for intermediate friction coefficients. The tail of the 
probability distributions depend also on friction. While 
P{AFn) decays exponentially for small frictions, it de- 
cays faster (slower) than exponential for intermediate 
(large) frictions (see the insets of Fig. 3]). 

Finally, we investigate the spatial distribution of the 
indeterminacy throughout the system. The aim is to find 
whether contacts that are located in a main force chain 
carry also larger force indeterminacies. In Fig. [5] (full 
circles) we plot the location of contacts that have larger 
Fn than twice of the average normal contact force (F„) 
(according to the original force network in the packing). 
We also plot the contacts with large force-indeterminacy 
AFn (open circles) above the average (AFn). This way 
approximately the same number of open and full cir- 
cles are plotted. It can be seen that contacts in force 
chains tend to have larger force indeterminacy in case 




FIG. 6: The correlation between the normal component of the 
original contact force Fn and the force indeterminacy A_F„ in 
terms of the friction coefficient ^. 



of intermediate friction [Fig. [Sfa)], but for ^ = 10 the 
two patterns become seemingly different [Fig. [S^b)]. In- 
deed, if we determine the correlation between lS.Fn and 
Fn over all contacts and plot it against the friction co- 
efficient (Fig. [B]), it reveals that the correlation vanishes 
for large frictions. Interestingly, the correlation exhibits 
again a nonmonotonic dependence on friction, where sig- 
nificant correlations are present for the intermediate fric- 
tion regime and weaker correlations outside. 

We note that local force indctcrminacies can be seen 
everywhere in our packings. This is in contrast to what 
has been reported in Ref. 0, where in case of large 
friction undetermined contacts formed localized clusters 
while contact forces in the rest of the packing became 



uniquely determined. Such localization is absent in the 
present study. We think that this difference originates 
from the different boundary conditions. In [3| the forces 
were kept fixed at the boundary which furthers the for- 
mation of a fully determined region. This is not the case 
here, where we prescribe only the global pressure. 

Conclusion - In this paper we presented the numer- 
ical results of the measurement of force indeterminacy 
in packings of frictional hard disks. We quantified the 
global force indeterminacy 77 of the packing with differ- 
ent methods and systematically studied the effect of in- 
terparticle friction coefRcient. 77 depends nonmonoton- 
ically on friction. We showed that the extent of force 
indeterminacy has an important influence on the me- 
chanical response properties of the material. The in- 
determinacy was also studied locally by measuring the 
interval of possible contact forces at individual contacts. 
We investigated the probability distribution of the in- 
tervals P(AF„) and the spatial distribution of the local 
indeterminacies. We found nonmonotonic friction depen- 
dence in the shape of P(AF„) and also in the correlation 
with contact forces. We observed significant correlation 
between the spatial pattern of the force-indeterminacy 
and force chains for intermediate frictions, however, the 
correlation disappeared in the large friction limit. 
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